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Wigner  symbols,  quantum  dynamics,  and  the  kicked  rotator 

J.  H.  Jensen 

Xiels  Bohr  Institute 
Blegdamsvej  17 
2100  Copenhagen  0,  Denmark 
and 

Department  of  Physics 
University  of  California 
Santa  Barbara,  CA  93106 

and 

Q.  Niu 

Department  of  Physics 
University  of  California 
Santa  Barbara,  CA  93106 

ABSTRACT:  Using  the  kicked  rotator  as  an  example,  we  show  how  to  apply  the 
formalism  of  Wigner  symbols  to  the  study  of  quantum  dynamics.  This  approach 
provides  a  simple  and  direct  method  of  comparing  quantum  and  classical  dynamics. 
We  investigate  in  detail  the  leading  quantum  effects  for  small  values  of  h  and  discuss 
the  time  scales  at  which  quantum  effects  appear. 


PAC'S  Indices:  0-3.65.-w,  0.3.65.Sq,  0.5.45.+b,  0.2.90.+p 


I.  INTRODUCTION 


Progress  in  understanding  nonin tegrable  classical  dynamical  systems1  has,  in  recent 
years,  stimulated  considerable  interest  in  the  corresponding  quantum  systems.  In  par¬ 
ticular,  much  work  has  been  devoted  to  investigating  relationships  between  the  classical 
and  quantum  dynamics  of  such  systems.  Important  subjects  of  investigation  include  the 
qualitative  differences  between  a  system's  classical  and  quantum  behaviors-,  the  time  scale 
at  which  quantum  effects  become  significant,  and  how  quantum  effects  can  be  calculated 
for  small  values  of  h? 

Insight  into  these  issues  may  be  found  by  using  a  formulation  of  quantum  mechan¬ 
ics  that  resembles,  as  much  as  possible,  the  usual  formulation  of  classical  mechanics.  A 
conventional  method  of  doing  this  is  to  construct  wave  packets  that  are  initially  localized 
around  a  particular  position  and  momentum.3  The  evolution  of  such  a  wave  packet  is  then 
compared  to  the  classical  trajectory  beginning  at  the  same  position  and  momentum.  While 
this  approach  has  proved  useful  for  many  applications,  we  advocate,  in  this  paper,  an  al¬ 
ternative  approach  that  we  believe  can  give  a  simpler  and  clearer  comparison  of  quantum 
and  classical  dynamics. 

The  central  idea  is  to  directly  calculate  the  time  dependence  of  quantum  operators, 
using  the  phase  space  representation  associated  with  the  works  of  Wigner  and  Weyl.*'* 
For  very  small  h,  quantum  operators  represented  in  this  manner  (sometimes  referred  to  as 
Wigner  symbols)  reduce  precisely  to  their  classical  counterparts  plus  corrections,  which  are 
typically  of  0{h“).  Wigner  symbols  thus  provide  a  quantum  generalization  for  any  classi¬ 
cal  quantity  that  has  an  associated  quantum  operator.  The  differences  between  quantum 
and  classical  dynamics  may  be  distinguished  by,  for  example,  comparing  the  Wigner  sym¬ 
bols  for  the  position  and  momentum  operators  with  the  classical  position  and  momentum 
trajectories. 
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The  use  of  Wigner  symbols  has  three  principle  advantages  over  the  use  of  wave  packets. 
First,  in  using  Wigner  symbols  the  only  parameter  added  to  the  classical  parameters  is  h. 
In  contrast,  with  wave  packets  one  must  also  specify,  at  least,  their  initial  widths.  Second, 
wave  packets  usually  spread  out  and  eventually  axe  no  longer  sharply  localized  around 
a  particular  position  and  momentum,  making  dubious  a  comparison  with  an  individual 
classical  trajectory.  With  Wigner  symbols,  this  spreading  problem  does  not  exist.  Third, 
by  employing  a  phase  space  representation  for  quantum  operators  the  difference  between 
quantum  and  classical  mechanics  may  be  separated  into  two  parts.  One  part  arises  from  the 
the  difference  in  the  phase  space  distribution  functions  allowed  for  quantum  and  classical 
systems  and  is  not  essentially  related  to  the  dynamics.  The  other  part,  that  which  is 
obtained  by  calculating  Wigner  symbols,  contains  all  the  dynamical  information.  Making 
this  separation  can  clarify  the  origin  of  quantum  effects. 

While  the  basic  formalism  of  Wigner  symbols  is  well-established,6  it  has  not,  to  our 
knowledge,  been  previously  applied  to  the  calculation  of  dynamical  properties  for  nonin- 
tegrable  quantum  systems.  In  this  paper,  we  demonstrate  how  to  do  this,  using  as  an 
example  the  kicked  rotator,  a  model  frequently  employed  in  the  investigation  of  noninte- 
grable  dynamics.'  This  model  consists  of  a  particle  that  is  periodically  “kicked”  and  is 
governed  by  the  Hamiltonian 

2  +x 

H  =  —  +  A' cos(x)  ^  6{t-n),  (1.1) 

n=-oc 

where  x  is  the  particle’s  position,  p  is  the  particle's  momentum,  and  K  is  the  kicking 
strength.6 

For  the  kicked  rotator,  we  show  explicitly  how  to  calculate  Wigner  symbols,  and  we 
study  in  detail  the  leading  deviation  for  small  h  of  symbols  from  their  classical  limits.  We 
find  that  the  asymptotic  time  dependence  of  the  leading  deviation  depends  only  on  general 
features  of  the  classical  trajectory.  We  also  discuss  how  to  estimate  the  time  scale  at  which 
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quantum  effects  begin  to  be  important  and  present  a  calculation  of  a  so-called  quantum 
crossover  time.  The  object  of  this  paper  is  to  demonstrate  the  utility  of  Wigner  symbols  for 
studying  quantum  dynamics,  as  well  as  to  give  new  results  for  the  kicked  rotator  relevant 
for  small  values  of  h. 

We  review,  in  Section  II.  the  basic  formalism  of  Wigner  symbols.  In  Section  III,  we 
give,  for  the  kicked  rotator,  a  practiced  method  for  calculating  Wigner  symbols  and  present 
typical  numerical  results,  which  we  compare  to  the  corresponding  classical  ones.  We  then, 
in  Section  IV,  consider  in  detail  the  leading  quantum  correction  to  a  classical  trajectory 
and  discuss  the  time  scales  at  which  quantum  effects  become  significant. 

II.  THE  WIGNER  SYMBOL  FORMALISM 


For  a  system  consisting  of  a  single  particle,  a  dynamical  quantity  A  =  F(x,p)  (e.g., 
the  particle's  energy)  can  classically  be  regarded  as  a  function,  A(xg,pg,t)  =  F(x(<),p(t)), 
of  the  particle's  initial  position  xq,  initial  momentum  po,  and  the  time  t.  Similarly,  a 
quantum  operator  .4  =  F(x,p)  can  be  regarded  as  a  function,  .4(xg,po,<)  =  F(x(t),p(t)), 
of  the  initial  position  operator  ig  =  x(0),  the  initial  momentum  operator  pg,  and  time. 
We  assume  that  A  is  Hermit ian,  that  the  function  F(x,p)  does  not  have  an  explicit  h 
dependence,  and  that  [x.p]  =  ih. 


Since  the  position  and  momentum  operators  do  not  commute,  the  function  A  is  well- 
defined  only  if  an  operator  ordering  is  specified.  Weyl  ordering  is  obtained  from  the 
prescription 


Mt)  =  J  <U*dn  Alro-Po,') J  ««(i(0)-..l+iTWi»-p.l,  (2.1) 

where  d  is  the  dimension  of  the  system.  Defined  in  this  manner,  we  refer  to  A  as  the 
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Wigner  symbol  for  A .®  To  find  A,  one  takes  the  Wigner  transform  of  A: 


A(xQ,p0,t)  =  J  da  e  *Qp°  (x0  +  a/2|  i(t)  |x0 -a/2),  (2.2) 

where  the  state  |i)  is  an  eigenstate  of  the  initial  position  operator.  The  Wigner  symbol  of 
any  Hermitian  operator  is  real. 


This  way  of  representing  a  quantum  operator  is  useful  because  the  Wigner  symbol 
is  a  natural  generalization  of  the  classical  A(io,po>0-  For  small  ft,  one  can 

show^ 


A{*0,P0J)  =  A(XQ,PQ,t)  +  0(ft2). 


(2.3) 


Furthermore,  the  average  of  .4  over  an  initial  state  |t/>o)  may  be  expressed  as 

(A(t))  =  J  dxdx'(ti>0\x)(x\A(t)\x')(x'\il>o) 

=  / ^'0^0--i(-ro>PO>Op(^0,Po)» 

where  p  is  the  Wigner  function  of  given  by4 

P(xO«Po)  =  J  e_^QP°(TO  +  a/2|0o)(^olIO-«/2). 


(2.4) 


(2.5) 


Equation  (2.4)  is  very  analogous  to  the  classical  expression  for  the  average  value  of  A  over 
a  distribution  of  initial  conditions. 


In  Eq.  (2.4),  the  ft-dependence  of  ,4’s  expectation  value  has  been  divided  into  two  parts, 
one  due  to  the  Wigner  function,  p ,  and  one  due  to  the  Wigner  symbol,  A.  The  Wigner 
function,  as  we  have  defined  it,  is  independent  of  time  and  has  no  connection  to  the  system’s 
Hamiltonian.  It  merely  represents  choice  of  initial  state.  Its  ft-dependence  arises  from  the 
basic  rules  of  quantum  mechanics  that  constrain  the  possible  forms  of  p.  For  example,  p 
cannot  be  choosen  to  be  a  delta  function  in  both  the  position  and  momentum  variables,  as 
this  would  violate  the  Heisenberg  uncertainty  relation.  The  Wigner  symbol,  on  the  other 
hand,  depends  essentially  on  the  Hamiltonian  and  contains  dynamical  information  intrinsic 
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to  the  system.  By  comparing  Wigner  symbols  directly  with  their  classical  limits,  one  may 
distinguish  quantum  effects  particular  to  a  system’s  dynamics  from  ones  due  solely  to  the 
general  restrictions  quantum  mechanics  places  on  phase  space  distributions. 

The  phase  space  formulation  of  quantum  mechanics  developed  here  is  distinct,  but 
complementary  to,  one  which  uses  time  dependent  Wigner  functions  as  the  central  dy¬ 
namical  objects.10  Our  approach  corresponds  to  the  Heisenberg  picture  of  quantum  me¬ 
chanics,  while  the  use  of  time  dependent  Wigner  functions  corresponds  to  the  Schrodinger 
picture.  In  previous  studies  of  quantum  dynamics,  the  Schrodinger  viewpoint  has  predom¬ 
inated.  In  the  following  sections,  we  show  how  practical  calculations  can  be  done  within 
the  Heisenberg  framework. 

III.  CALCULATING  WIGNER  SYMBOLS 
FOR  THE  KICKED  ROTATOR 


The  essential  features  of  the  kicked  rotator's  dynamics  are  contained  in  a  map  that 
gives  the  position  and  momentum  just  before  a  kick  in  terms  of  the  position  and  momen¬ 
tum  before  the  previous  kick.  This  map  is  found  by  integrating  the  equations  of  motion 
corresponding  to  the  Hamiltonian  (1.1)  over  one  kicking  period,  leading  to  the  so-called 
standard  map1 1 


•^’n+l  —  xn  "b  Pn+1 
p»i  +  l  =  Pn  +  sin(j„), 

with  in  being  the  position  and  pn  the  momentum  before  the  nth  kick. 


(3-1) 


The  dependence  of  a  dynamical  quantity,  .4.;J  =  F(xn,Pn)i  on  the  initial  position  and 
momentum  can  be  obtained  from  the  recurrence  relation 


.4;i+[(.r0.po)  =  --M-fO  +  P0  +  A'sin(xo),po  +  Ksin(i0)),  (3.2) 


which  follows  from  (3.1)  and  the  fact  that  an  (n  +  l)-step  evolution  from  (xfbPo)  *s  the 
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same  as  an  n-step  evolution  from  (j^pi).  Moreover,  by  applying  (3.2)  to  the  cases  with 
-4n  =  xn  and  An  =  pn,  (3.1)  can  be  derived  from  (3.2).  Thus  (3.1)  and  (3.2)  provide 
equivalent  representations  of  the  kicked  rotator’s  classical  dynamics. 

Quantum  mechanically,  the  evolution  of  an  operator  over  one  kicking  period  is  given 

byl0.12 

An+I  =  U'AnU,  (3.3) 

where  U  is  an  evolution  operator  defined 

U  =  e~k(rt/2)e-*K  cos(io).  (3.4) 

Using  Eq.  (2.2)  to  express  (3.3)  in  terms  of  Wigner  symbols,  one  finds,  after  some 
algebra, 

-t-CC 

.4n+1(xo,po)  =  ^  Jl[2Ks\n(xQ)/h]An(xo  +  pQ  +  hl/2,po  +  hl/2)>  (3.5) 

/=-* 

where  J/( x)  is  a  Bessel  function.  The  recurrence  relation  (3.5)  is  a  quantum  generalization 
of  (3.2),  reducing  to  (3.2)  in  the  k  — *  0  limit.  For  large  /,  the  Bessel  functions  decrease 
rapidly  and  only  about  4 K/h  terms  contribute  substantially  the  sum  in  (3.5). 

The  numerical  evaluation  of  (3.5)  is  especially  simple  for  quantities  with  the  symmetry 
An(*0'P0)  =  -4nUo  +  2?r,po)  =  .4n(x0,po  4-  2?r),  (3.6) 

which  is  preserved  by  the  dynamics.  This  symmetry  can  be  exploited  by  choosing  h  equal 
to  a  "resonant”  value  of  4 Ka/b,  where  a  and  b  axe  positive  integers.'  12  Any  other  value  of 
h  may  be  approximated  to  arbitrary  accuracy  by  using  sufficiently  large  a  and  6. 13 

For  these  resonant  values  of  h,  (3.5)  becomes 

6-1 

>Wl(-r0iP0)  =  £^0Mn(*0  +po  +  ^/2,P0  +  hi/ 2),  (3.7) 

1=0 
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with 

\  1  V"  [ 2kIju  K.  2irm.  ,  2k m  ' 

4?K*0)  =  £  2^  cos  ~ ~  +  y  («*(*<>  +  -£- )  -  cos(x0 - — ))  ,  (3.8) 

m=0 

which  now  contains  a  finite  number  of  terms.  Further  simplification  occurs  if  pg  is  ap¬ 
proximated  by  a  number  2 Kj/b.  where  j  is  am  integer.  The  recurrence  relation  (3.7)  is 
then  equivalent  to  a  map  of  a  6x6  matrix  and  is  straightforward  to  evaluate  numerically. 
The  symmetry  (3.6)  applies  in  particular  to  the  momentum  step,  pn+l(*0>Po)  —  pn(xo,Po), 
from  which  the  momentum  symbol  can  be  obtained.  In  fact,  the  Wigner  symbol  for  any 
operator  that  is  a  polynomial  in  x  and  p  can  be  found  by  iterating  6x6  matrices. 

In  Fig.  1,  we  show  the  time  dependence  of  a  typical  momentum  symbol  for  two  values 
of  h ,  as  well  as  the  corresponding  classical  trajectory  ( K  =  2.5,  xq  =  1.0,  pg  %  0.94).  The 
symbols  were  calculated  by  using  (3.7)  with  6  =  67.  The  classical  trajectory,  Fig.  la,  is 
chaotic.  For  small  h.  the  momentum  symbol  also  appears  irregular,  as  in  Fig.  lb,  but  as 
h  is  increased  it  becomes  more  regular  as  suggested  by  Fig.  lc. 

A  Fourier  transform  in  time  of  the  momentum  symbol  consists  of  sharp  peaks,  indicat¬ 
ing  quasiperiodic  dynamical  behavior.  For  small  h.  many  peaks  have  a  substantial  weight, 
while  for  large  h  only  a  few  are  significant.  In  the  h  — *  0  limit,  the  number  of  important 
Fourier  peaks  may  diverge,  leading  to  chaotic  classical  motion.  The  suppression  of  the 
chaotic  behavior  of  the  kicked  rotator  by  quantum  effects  can  be  understood  in  terms  of  a 
localization  argument  due  to  Grempel.  et  al.14 

IV.  QUANTUM  CORRECTIONS  FOR  SMALL  h 

We  now  consider  the  leading  deviation  of  a  Wigner  symbol  for  the  kicked  rotator  from 
its  classical  limit.  These  deviations,  generally  of  0(h 2),  represent  the  first  quantum  cor¬ 
rections  to  the  classical  motion.  Here  we  restrict  ourselves  to  the  position  and  momentum 
symbols,  although  the  method  we  use  may  be  extended  to  other  Wigner  symbols. 
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The  position  and  momentum  symbols  have  the  expansions 


*n(-ro-Po)  =  -rnUo-Po)  +  h2an{x  o,po)  +  0(&4), 


(4.1) 

P»il*0-P0)  =  P/iUu-Po)  +  fi2Ai(*0,po)  + 

where  an  and  0n  axe  the  coefficients  of  the  leading  quantum  corrections.  These  coefficients 
may  be  obtained  by  iterating  the  map 


—  +  ,-^n+l 

3n+l  =  3n  -  vW(xn)an  +  Wn, 
with  the  initial  condition  oq  =  Jq  =  0-  (4.2), 


(4.2) 


I\  1  X 

Wn  =  j€ik€jixn,ij{-sm(xn)xn<iei  +  -cos(j  n)xn>kxnj},  (4.3) 

where  we  use  the  notations 

/j(xo,W) -2figjS2l. 

/.»<* O.PO)  (4'4) 

/.|>(J0»P0)  =[/,i(J'0-Po)j.;- 

for  any  function  /(xo-Po)-  In  (4.3),  the  indices  z,^.  A:.  /  are  summed  from  1  to  2,  and  e(j  is  an 
antisymmetric  matrix  with  ep  =  «22  =  0  and  =  — *21  =  T  Wn  can  be  calculated  from 
the  classical  map  (3.1).  While  (4.2)  may  be  derived  from  (3.3).  a  more  direct  derivation  is 
given  in  the  Appendix. 


For  short  times,  the  quantum  corrections  obtained  from  (4.2)  are  small,  and  the  quan¬ 
tum  and  classical  trajectories  nearly  coincide.  However,  as  time  increases,  the  quantum 
corrections  tend  to  grow  and  eventually  become  important  (<ee  Fig.  2).  For  small  h,  it  is 
the  asymptotic  behavior  for  large  time  of  the  corrections  which  is  of  interest.13  Numerical 
study  of  (4.2)  indicates  the  asymptotic  form  of  the  corrections  depends  only  on  the  general 
character  of  the  classical  trajectory.  Specifically,  we  find  for  a  generic  quasiperiodic  clas¬ 
sical  trajectory  that,  on  the  average,  |an|  and  |/in|  increase  proportionally  to  n3,  while  for 
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a  chaotic  classical  trajectory  they  increase  proportionally  to  e3<7T\  where  a  is  the  largest 
Liapunov  exponent  for  the  classical  trajectory.  Equations  (4.1)  then  suggest  that  the  time 
at  which  the  Wigner  symbols  begin  to  deviate  significantly  from  the  corresponding  classical 
trajectories  scales  as  h~ 2/3,  for  quasiperiodic  trajectories,  and  as  —  (2/3<r)  \n(h)  for  chaotic 
trajectories.16  Clearly,  the  quantum  corrections  to  a  chaotic  classical  trajectory  increase 
much  more  rapidly  than  those  for  a  quasiperiodic  trajectory. 

Figure  3  gives  a  comparison  of  the  quantum  corrections  computed  from  (4.2)  with  exact 
results  obtained  from  Eq.  (3.7)  for  two  trajectories  with  K  —  1.5  and  h  =  47r/250  «  0.05. 
We  define  the  quantum  deviation  An  =  (pn  -  Pn)A2-  For  small  n,  An  is  approximately 
equal  to  dn,  but  as  n  increases  the  deviation  grows  and  can  no  longer  be  obtained  from 
(4.2).  The  deviation  for  the  chaotic  trajectory  (squares)  initially  grows  much  faster  than 
the  deviation  for  the  quasiperiodic  trajectory  (circles). 


As  an  illustration  of  how  to  apply  these  results,  we  consider  a  wave  packet  that  initially 
has  the  Gaussian  form 

=  (A)  '  exp  -  A(*  -  i0)2),  (4.5) 

where  A  is  a  real  parameter.  From  Eq.  (2.5),  one  finds  the  corresponding  Wigner  function 


to  be 


PU.P)  =  A  exp  (A(*  -  x„)2  -  ^P-po)2)- 


The  average  position  of  this  wave  packet  as  it  evolves  in  time  is,  as  follows  from  (2.4), 


{x(t))  =  J  dx'dp'  £(x',p',t)p(x'.p). 


We  now  ask  how  this  average  position  differs,  for  small  h,  from  the  average  position 
found  from  the  classical  evolution  of  a  phase  space  distribution  that  is  initially  set  equal 
to  the  Wigner  function  (4.6).  In  other  words,  we  compare  the  wave  packet’s  average 
position  with  that  obtained  from  Eq.  (4.7)  when  x  is  replaced  by  the  classical  position 
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X. 


Using  (4.1)  one  finds  that  the  difference  between  these  quantum  and  classical  average 
positions  is  simply  7i2an(xo,po)  plus  terms  of  higher  order  in  h.  Hence,  the  time  at  which 
the  average  position  of  the  Gaussian  wave  packet  begins  to  deviate  substantially  from  the 
classical  average  is  the  same  as  the  time  at  which  xn(io,po)  begins  to  deviate  substantially 
from  x fj ( x 0 ,  po ) •  We  note  that  this  time  cannot  be  obtained  from  the  usual  semiclassical 
treatment  of  Gaussian  wave  packets,  since  the  average  position  in  this  approximation  is 
just  the  classical  xn(xo-Po)-l‘ 

As  a  further  example,  we  apply  Wigner  symbols  to  estimate  the  so-called  quantum 
crossover  time  for  the  kicked  rotator.  The  crossover  time  may  be  defined,  somev  ,iat  loosely, 
as  the  time  at  which  there  appears  a  significant  qualitative  difference  between  the  kicked 
rotator’s  classical  and  quantum  behavior.  For  values  of  K  greater  than  Kc  =  0.9716..., 
this  is  conventionally  taken  as  the  time  at  which  the  average  momentum  of  a  quantum 
wave  packet  ceases  to  exhibit  the  diffusive  behavior  observed  classically.7  (If  I\  is  less  than 
A'c,  classical  KAM  trajectories  bound  the  momentum  and  this  criterion  does  not  apply.11) 

The  classical  diffusion  rate  for  K  >I\ ' c  is  determined  essentially  by  the  dynamics  in  the 
vicinity  of  the  cantorus  remnant  of  the  last  bounding  KAM  trajectory  to  disappear  (which 
it  does  at  K  =  Kc ) •  Since  this  last  KAM  trajectory  plays  a  central  role  in  determining 
the  kicked  rotator's  qualitative  behavior,  it  is  reasonable  to  postulate  that  the  time  at 
which  quantum  corrections  to  this  trajectory  become  significant  should  correspond  to  the 
quantum  crossover  time.  Numerical  calculations  obtained  from  Eq.  (4.2)  for  the  leading 
quantum  corrections  to  the  position  and  momentum  symbols  indicate  that  these  grow 
proportionally  to  nu.  where  v  =  6.2  ±0.1  (see  Fig.  4).l!*  Equations  (4.1)  then  imply  a 
crossover  time  proportional  to  which  is  consistent  with  both  numerical  simulations 

and  with  the  result  of  Fishman,  et  al., 19  who  use  a  renormalization  argument  to  find  a 
crossover  time  scaling  as  h  ^2/*/  ),  with  v'  «  6.1. 
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V.  CONCLUDING  REMARKS 


Both  classically  and  quantum  mechanically,  there  are  two  complementary  perspectives 
from  which  one  may  discuss  dynamics.  The  usual  classical  approach  is  to  calculate  trajec¬ 
tories  for  individual  particles.  Alternatively,  one  may  study  the  evolution  of  a  distribution 
function  of  phase  space  points,  by  solving,  for  example,  Liouville’s  equation.  In  quantum 
mechanics,  the  viewpoint  analogous  to  the  classical  trajectory  approach  is  the  Heisenberg 
picture,  while  the  Schrodinger  picture  is  analogous  to  the  classical  distribution  function 
approach. 

Since  many  of  the  important  concepts  pertaining  to  nonin tegrable  classical  dynamical 
systems  (e.g.,  the  distinction  between  regular  and  chaotic  dynamics)  are  most  easily  un¬ 
derstood  from  the  trajectory  viewpoint,  it  seems  natural  to  employ  the  Heisenberg  picture 
in  comparing  quantum  and  classical  dynamics.  Nevertheless,  this  is  rarely  done,  perhaps 
because  of  a  lack  of  practical  techniques  for  applying  the  Heisenberg  approach.  In  this 
paper,  we  have  sought  to  develop  such  techniques,  using  the  model  of  the  kicked  rotator. 

Our  basic  method  is  to  calculate  directly  the  time  dependence  of  the  Wigner  symbols 
corresponding  to  quantum  observables.  When  quantum  mechanics  is  formulated  in  this 
way,  its  relationship  to  classical  mechanics  is  very  clear.  Each  Wigner  symbol  corresponds 
to  a  unique  classical  quantity,  allowing  the  comparison  between  quantum  and  classical 
dynamics  to  be  made  simply  and  unambiguously.  For  example,  the  Wigner  symbol  cor¬ 
responding  to  the  derivative  of  a  classical  trajectory  with  respect  to  its  initial  position, 
dx{t)/dx q.  which  is  relevant  to  the  calculation  of  Liapunov  exponents,  is  simply  the  deriva¬ 
tive  of  the  position  symbol,  dx(t)/diQ.  How  to  make  such  a  connection  in  the  language  of 
Schrodinger  wave  functions  is  less  obvious. 

In  addition  to  giving  a  practical  algorithm  for  calculating  Wigner  symbols  for  the  kicked 
rotator,  we  have  also  studied  the  behavior  of  the  position  and  momentum  symbols  for  small 


12 


values  of  h.  We  find  that  the  asymptotic  behavior  of  the  leading  quantum  corrections  for 
these  Wigner  symbols  depends  only  on  the  general  character  of  the  corresponding  classical 
trajectory,  and  we  have  shown  that  these  corrections  can  be  used  to  find  the  time  at  which 
quantum  effects  become  significant.  In  particular,  the  quantum  crossover  time  for  the 
kicked  rotator  with  K>KC  can  be  estimated  from  the  deviations  to  the  last  bounding 
I<AM  trajectory. 

While  in  this  paper  we  have  r<  icted  ourselves  to  the  kicked  rotator,  the  Wigner 
symbol  approach  can  be  applied  to  many  other  systems.  This  method  is  especially  useful 
when  a  comparison  of  classical  and  quantum  dynamics  is  desired.  We  anticipate  applica¬ 
tions  both  to  questions  concerning  the  qualitative  dynamical  behavior  of  simple  quantum 
systems  and  to  the  development  of  practical  quantitative  calculational  techniques. 
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APPENDIX 

Here  we  present  a  derivation  of  map  (4.2),  which  gives  the  leading  quantum  corrections 
for  the  position  and  momentum  symbols.  We  begin  with  some  general  considerations 
pertaining  to  the  h  expansions  of  Wigner  symbols.  In  developing  such  expansions,  it  is 
convenient  to  have  a  relation  between  the  Wigner  symbols  for  two  individual  operators,  A 
and  5,  and  the  Wigner  symbol  for  their  product,  AB  =  C .  This  is  given  by  the  so-called 
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‘‘♦-product:” 

/>/  ,  -  \ih  (  d  d  d  d  \  - 

C(x„,P0)  =.4(x0.M)exp[7(^^  -  B(x„,p0)  (  4  „ 

=  A(xQ.pQ)  *  B(sq,po ), 

with  the  arrows  indicating  on  which  sides  the  derivatives  act.6,20  The  *-product,  like  the 
operator  product,  is  associative  but  not  commutative.  Note  that  the  *-product  reduces  to 
an  ordinary  product  for  h  —  0. 


Expanding  the  exponential  in  Eq.  (A.l),  we  obtain  the  h  expansion 


00  i  /  ih\n 

A(*0<P0)  *  B(T0,b0 )  =  E  ?  {**>.. 

«— ft  '  ' 


(A.2) 


where 


rS(  }  (^)a(n-^xoa-po^x0a(«-m)po‘ 

{A,B}  1  is  simply  the  conventional  Poisson  bracket.  Evaluating  a  *-product  of  two  Wigner 


(A.  3) 


symbols  to  order  n  in  h  requires  derivatives  up  to  order  n. 


We  now  consider  the  h  expansion  for  the  Wigner  symbol  of  some  analytic  function  of 
an  operator.  F(A).  The  Wigner  symbol  for  F(A)  may  be  written  F[*A(xq, po)],  where 
F[*A(xo,po)]  is  defined  in  terms  of  F’s  Taylor  expansion  and  the  rules 

*A(j0-Po)  — A(xO'Po)' 

[*A(x0,po)]2  =A(T0- Po)  *  A(.ro.po),  (A.4) 

[*A(x0.po)]3  =.4(.r0,po)  *  .4(xo«Po)  *  A(x0,po),  etc. 

We  may  expand  F[*.4(xq.po)]  as 

F[*.4(.r0.po)]  =F[.4(.r0,po)  -M*A(x0,Po)  “  A(x0,po))] 

=  V  — yF(,,)[.4(xo,po)][*A(x0,po)  -  A(x0,po)]n, 

n=0  n’ 

where  F^n\A)  is  the  nth  derivative  of  F(A).  The  n  =  i  term  varnishes  identically.  Using 
(A. 2),  each  term  of  the  right  hand  side  of  (A. 5)  can  be  further  expanded  in  powers  of  h. 
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To  o(n2),  one  finds  that  only  the  first  three  nontrivial  terms  (i.e.,  n  =  0,2,3)  contribute, 
with  the  higher  order  terms  being  of  Oih*)  or  smaller.  Evaluating  these  gives 
F[*i(x0,p0)]  =F[.4(.r0.p0)] 

~  y«z/fcej/^.Ijl-r0!P0){^(2)[^(x0,P0)]^)H(I0»P0)  (A6) 

+  ^F(3)[ii.r0.po)]i.4.(xo,po)i/(xo»Po)}  +  0(h4), 

where  we  have  used  the  notations  of  (4.4). 

The  dynamics  of  the  quantized  kicked  rotator  follows  from  the  operator  map 

*n+l  =  in  +Pn+1  . 

(A.  7) 

Pn+L  -  Pn  +  K  sm(ln), 


which  is  the  quantized  form  of  (3.1 )  and  is  equivalent  to  (3.3).  Taking  the  Wigner  transform 


of  this  map  gives 


J’n+l  —  r'i  4-  pnjf  i 


(AS) 

Pn+l  =  Pn  +  A'  sin(*i„). 

Applying  (A. 6)  and  (4.1)  to  (A.S)  and  extracting  the  terms  of  0(7i2),  leads  directly  to 
(4.2). 
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FIGURE  CAPTIONS 


1)  A  classical  momentum  trajectory  and  its  corresponding  momentum  symbol  for  K  = 

2.5,  iq  =  1.0,  and  p q  =  20:r/67:  (a)  is  the  classical  trajectory,  (b)  is  the  momentum  symbol 
with  h  =  8t/67  ss  0.38.  and  ( c )  is  the  symbol  with  h  =  2127t/67  as  9.9.  As  h  is  increased 
the  momentum  symbol's  time  dependence  evolves  from  chaotic,  in  the  classical  limit,  to 
simple  and  quasiperiodic,  for  large  values  of  h. 

2)  Momentum  symbol  (h  =  4t/250;  solid  line)  with  the  corresponding  classical  trajec¬ 
tory  (dashed  line)  for  K  =  1.5,  iq  =  3.0,  and  pq  =  n.  The  momentum  symbol  follows  the 
classical  trajectory  for  several  iterations,  but  then  begins  to  deviate. 

3)  Logarithm  of  the  deviation  An  s  (pn  -pn)/ft2  for  a  quasiperiodic  trajectory  ( K  = 

1.5.  xq  =  3.0.  and  pq  =  r:  circles)  and  a  chaotic  trajectory  ( K  —  1.5,  io  =  2.0.  and  po  =  t; 
squares).  The  solid  lines  are  exact  calculations  with  h  —  47r/250,  rind  the  dashed  lines  are 
the  approximations  obtained  by  using  the  h  expansion  coefficient  /3n. 

4)  Logarithm  of  the  quantum  correction  an  for  the  last  bounding  KAM  trajectory 
(K  =  0.971G35 . . ..  £Q  =  7r,  and  po  =  3.73799S . . .).  Only  points  from  the  curve’s  upper 
envelope  are  plotted  (solid  line).  The  dashed  line  is  a  least  squares  fit  of  slope  6.2,  showing 
that  the  correction  increases  approximately  as  n6  ‘2. 
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